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Abstract 

From a common expression for the poloidal electrostatic field of a tokamak, in the limit of large 
aspect ratio and concentric circular flux surfaces, one may determine the associated potential. This 
potential satisfies Poisson's equation, which reduces to Laplace's equation when the medium has 
vanishing charge density, in axial geometry but not toroidal geometry. A simple transformation 
takes the potential over to the correct harmonic form for tokamak coordinates, and the resulting 
electrostatic field is calculated. From the radial field one may estimate the supporting charge 
density on the boundary, and from the poloidal field one may determine a prediction for the radial 
dependence of the electron temperature, which does not compare well with a rough estimate of the 
profile often seen in a tokamak. 
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I. INTRODUCTION 



From a common expression for the poloidal electrostatic field of a tokamak, in the limit 
of large aspect ratio and concentric circular flux surfaces, one may determine the associated 
potential. This potential satisfies Poisson's equation, which reduces to Laplace's equation 
when the medium has vanishing charge density, in axial geometry but not toroidal geometry. 
A simple transformation takes the potential over to the correct harmonic form for tokamak 
coordinates, and the resulting electrostatic field is calculated. From the radial field one may 
estimate the supporting charge density on the boundary, and from the poloidal field one 
may determine a prediction for the radial dependence of the electron temperature, which 
does not compare well with a rough estimate of the profile often seen in a tokamak. 



II. DERIVATION 



We start with the expression given by Wesson |l( f° r the poloidal electrostatic field at 
equilibrium 



(E 4> B^/B e }B 2 E^Btf, , ( (l/B e }B 2 I x 

Ee= WmBe -^ + RB " Pm {(BVB e }B e -B-) ' (1) 

where (A) is the flux surface average of A. In the limit of large aspect ratio e < 1 and 

concentric circular flux surfaces Aqs(t) = 0, and assuming toroidal symmetry d/d<f) = 0, 

one may write Eg = ^^(l + e cos 6) — Ejp / (1 + e cos 6>), which reduces to {2] 

E e = Ef\±ecos9 -£ 2 cos2fl)/2 , (2) 

upon application of Stokes' theorem to Faraday's law [3]. The associated potential may be 
found from Eg = —d^/rdO, which evaluates to $ axi = $o(8£ sin # — e 2 sin26 l )/4. Note that 
this $0 is not the $° of the GT-FRC model 0, [5] but is a unit bearing constant which sets 



the scale. The reason for the subscript "axi" for "axia! 



becomes apparent when one notices 



that the potential is of the correct harmonic form [(j, [7] for an unweighted circle representing 
a cylindrical flux surface in (r, 9, z) coordinates with axial symmetry, as is easily verified 
in (Z = ~rsin6,R = Rq + rcos#, z) coordinates via application of the axial Laplacian 
V 2 axi = d 2 /dZ 2 +d 2 /dR 2 to $ axi = %Z[(R-R ) /2R 2 -2/ R ]. In order to achieve the correct 
harmonic form for tokamak geometry, the potential must satisfy the toroidal Laplacian, 
which in cylindrical coordinates (3], a [9] is given by V 2 or . = V 2 axi + d/RdR (note that the 



expression for A* given by Hopcraft in Dendy's textbook lOj is not the toroidal Laplacian 
and differs by the sign of the additional term), and the form of the additional geometric term 
hints at the solution. Direct integration yields the toroidal potential <& tor = & axi (R — > lni?), 
from which E z = -$ [(lni? - Rq)/2B% - 2/R ] and E R = -$> Z/2RR 2 , noting that the 
introduction of the logarithm breaks the usually obvious relation between the symbol for 
the magnitude of a quantity and the units associated with that quantity — carefully pulling 
the units into the leading coefficients of expressions ensures that they are respected. From 
these, we get the poloidal field 
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coordinates, we have $ tor = $o r sm9[— ln(i? + ^cos^) + 5i?o]/2-Ro; anc l 
$o f - r, , „ „s r 2 sin 2 
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As this is an electrostatic field supported by a neutral medium (ie one for which the net 
charge on a differential volume element vanishes), Maxwell's equations V x E = and 
V • E = are satisfied. 

As electrostatic fields necessarily require a supporting charge density 
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13. 



one may very well ask where these charges are. As the solution to Laplace's equation is 
uniquely determined by the boundary condition, we might as well put them on the boundary 
of the region under consideration, which in this case is the R/Ro weighted circle representing 
our outermost flux surface at normalized minor radius r = 1 upon collapse of the toroidal 
dimension, giving us an inverse Dirichlet problem, which is usually defined as solving for 
the potential given the boundary condition. As the current in a tokamak is observed in the 
toroidal direction, we suppose the plasma to behave as a dielectric in the plane orthogonal 
to the current, and take the polarization to be proportional to the electrostatic field, P oc E, 
an admittedly crude approximation but one which satisfies the requirement of no volumetric 
bound charge density, p b = —V ■ P = 0. Were the plasma conductive in the (Z, R) = (r, 9) 



plane, no electrostatic field could be supported Q|. To find the surface bound charge density, 
we take — P • r, noting that quantity represents the charge density on the weighted circle, 
so to get the physical charge density on a poloidal circle of the toroidal flux surface we take 
a b (TbRo/R- We note that usually it is an applied electric field which drives the polarization 
of dielectrics, but for a magnetized plasma other effects may support the polarization. 

From the expression for the poloidal field, one may determine a prediction for the as- 
sociated electron temperature profile by solving the flux surface Fourier moments of the 
poloidal equation of motion 5|. With electron density n e = n° e {l + n^cos^ + n^sin^) and 
temperature kT e = T e , and the unity, cosine, and sine moments defined by the expressions 
(A)u,c,s = § d0{l,cos0,sin0}(l + ecos9)A/2ir, we solve (T e dn e /rd6 + n e eE e )u,c,s = for 
(the pure numbers) n c e ,n% and Ef \ given by 



(7) 



where F = 5i?o — hi-Ro, i n terms of the functions G^tj) = [(6-F — l)e 2 — 8F] and 

G 1 (r) = y/(F + 2)[{16F[4F - (1 + 3F)e*]} + 9e 4 ] , (8) 

from which we obtain $ = T<$VQRoT e /eGi. Note that the quantities appearing in, between, 
and following Equations ([7]) and (IE]) carry no units, as the units cancel out of the system 
of equations solved. The factor of Ro appearing in the scale factor of $o has a physical 
interpretation — for an equivalent supporting charge density and minor radius, as one scales 
the system by Ro, the magnitude of the charge distribution on the weighted circle increases 
as the major circumference of the torus, thus increasing the magnitude of the potential by a 
numerical factor of Ro. As $0 is constant V r, the explicit r = eRo dependence of G\ must 
cancel the implicit r dependence of T e (r), allowing us to write T e {r) /G\{r) = 5 ? e (0)/Gi(0), 
thus Gi determines the r dependence of the electron temperature as a prediction of the 
hypothesis that the electrostatic field in a tokamak is of a form similar to that of Equation ([2]). 
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III. PRESENTATION 



The following presents the results of a numerical evaluation of the expressions in the 
preceding section. The scale is set by Ro = 2 and T e (0)/e = 1, and r is the normalized 
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minor radius. Some useful geometry is presented in the Appendix. The toroidal electrostatic 
potential, &tor, is shown in Figure 1 — note the correct harmonic form of a stretched circular 
membrane (eg a drumhead put on a hoop given by the boundary condition), consistent with 



the physics of potential theory 



7| , and the obvious lack of bumps or poles, consistent with 
a neutral medium (the term "pole" actually means something physical SB , ll5) , namely the 



effect on the potential of an isolated flux-source which ultimately is quantized in units of 
the charge). We stress that these are multiple views of the same single object. In Figure 2 
we display the magnitude of the associated electrostatic field, and in Figure 3 we display 
the electrostatic field and the associated supporting charge density. 

The normalized electron temperature profile in a tokamak often compares favorably (from 
unpublished results of the analysis behind Reference 4J) with the formula T e (r)/T e (0) « 
(1 — r 4 ) 2 , which misses the effect of the pedestal but works surprisingly well for r < .8, 
and so that is the assumed profile to which we will make our comparison in Figure 4. The 
normalized predicted profile for the temperature is G\(r) / G\(0) and fails to fall sufficiently 
quickly to match the shape of the expected temperature profile. Thus, while comparison with 
our approximate formula hardly qualifies as a true experimental test of the hypothesis [if]], 
the existence of an electrostatic field of form similar to that of Equation (j2J) is not supported 
by this analysis. 



IV. CONCLUSION 



By starting with the neoclassical prediction for the poloidal electrostatic field in a toka- 
mak, we have determined the corresponding electrostatic field appropriate for toroidal de- 
vices, as well as the supporting charge density on the boundary. Examination of the solution 
to the electron poloidal equation of motion reveals a prediction for the temperature profile, 
which does not compare well to that expected to be measured in a tokamak. 



APPENDIX: SOME USEFUL GEOMETRY 



There are (at least) three useful sets of coordinates to describe a toroidal magnetic 
confinement device with concentric circular flux surfaces, namely (Z,R,<f>), (r,9,<f>), and 
(r, _L, ||), and in the infinite aspect ratio limit (Rq — > oo) we have the axial coordinates 
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(Z,R,z) and (r,0,z). (Note that the term "toroidal coordinates" means something very 
different to a mathematician than those commonly applied to a tokamak, which we call 
"tokamak coordinates.") The relationship between vectors in the tokamak coordinates 
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where || = B/£> = b = (0,bg, b^). Various operators in tokamak coordinates are given by 
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and some useful relationships are 
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(A.2) 
(A.3) 

(A.4) 
(A.5) 



Z = -r sin 9 , OZ/Or = Z/r , dZ/d9 = -(R - R ) , 

R - R = r cos 9 , dR/dr = (R - R )/r , dR/d9 = Z , 



(A.6) 



and for the logarithm, we have 9 In R/OR{m) = 1/R(m), which shows that the lni? carries 
no units. 
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Electrostatic Potential 
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FIG. 1: (Color online.) The toroidal electrostatic potential <&tor seen from multiple views. 



Electrostatic Field Magnitude 
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FIG. 2: (Color online.) The magnitude of the associated electrostatic field seen from multiple 
views. 
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Electrostatic Field and Supporting Charge Density 
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FIG. 3: (Color online.) The electrostatic field and associated supporting charge density seen from 
multiple views. 



1 

0.9h 



fi°-8r 
f-0.7h 



0.6 
0.5 



.2 

03 

"8 °- 4 

| 0.3 
I 0.2 



0.1 









Prediction vs Approximation 
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FIG. 4: Comparison of the neoclassical prediction for the normalized electron temperature profile, 



G\, versus an approximation to that found in tokamak experiments, (1 
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